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Abstract. We establish an exact formula for the average number of edges appearing 
on the boundary of the global convex hull of n independent Brownian paths in the 
plane. This requires the introduction of a counting criterion which amounts to "cutting 
off" edges that are, in a specific sense, small. The main argument consists in a mapping 
between planar Brownian convex hulls and configurations of constrained, independent 
" linear Brownian motions. This new formula is confirmed by retrieving an existing 

exact result on the average perimeter of the boundary of Brownian convex hulls in the 
plane. 
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^ ■ Introduction 
Ov 

\ Consider n independent Brownian paths in the plane, all of the same given duration 

and starting from the same point, and let C n denote their global convex hull, that is 
the smallest convex set containing all n paths (figured]). As an example of randomly- 
produced convex shapes based on stochastic processes, C n (particularly when n = 1) 
has received continuous attention over the past decades, in the fields of probability 
theory [THT2] and stochastic geometry [T3l[Tl] . with applications ranging from image 
processing and pattern recognition [15] to statistical physics [Ej, ecology p]3,[TB] and 
biophysics [19]. (We focus on Brownian convex hulls, given the fundamental position 
that Brownian motion has been occupying in the mathematical, physical and biological 
sciences for more than a century [20H26] ; the literature on random polytopes is vast, 
and it is not our aim here to present even just an overview of it: we refer the reader 
to [2TH3T] for surveys.) 

Applications along with theoretical considerations lead to common, typical 
questions about random polytopes, pertaining for instance to the average values of 
some elementary geometrical quantities such as: perimeter, enclosed area or volume, 
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Figure 1: A realization of dCg (solid black line), the boundary of the convex hull of 
n = 9 independent planar Brownian paths, all starting from the origin and having the 
same duration. (The dashed lines are canonical x and y-axes.) Each planar Brownian 
path is defined as having for coordinates two independent linear Brownian motions. 

number of edges or sides. In the case of the global convex hull of n planar Brownian 
paths, exact formulae for the average perimeter and enclosed area exist in the literature: 
see [HIEIES] for n — 1 and [HUB] for n > 2. The n > 2 - case led to the design of a 
general method for the computation of average perimeters and areas, shedding light 
on the links between random convex hulls and extreme- value statistics [16]. However, 
somewhat finer aspects concerning the shape of these hulls do not seem amenable to 
this method. For instance, the intuitive conjecture that, for large n but a fixed, finite 
common duration, C n tends in some sense to a circular shape has been proven recently 
using other techniques [33]. We aim to pursue here the study of the shape of C n by 
examining the set of line segments which compose its boundary dC n : indeed, if Brownian 
convex hulls do not have any isolated extreme point [6], they consist almost surely of 
countably many line segments [SI El IS]- This suggests that, provided a regularization 
is introduced to "cut off" some of the "smaller" edges, one can compute the average 
number of edges on the convex hull of planar Brownian paths. 

In the mathematical literature, exact formulae for average numbers of edges, sides 
or vertices are to be found in contexts where the random convex hull is that of a discrete 
set of independent random points [31HH] • But similar results remain rare in the case of 
dependent points correlated by their being the successive positions of a random walk. 
Our starting point is one of these results: a combinatorial lemma by G. Baxter [12|B3] 
leading to an exact formula for the average number (V(N)) of edges on the convex hull 
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of an iV-step random walk in the plane: 

N 1 

( V (N))=2J2t- (1) 

fc=i 

In this paper, we first show that, with a proper regularization, Baxter's formula 
extends to the convex hull of a single planar Brownian path precisely as one would 
expect by taking the continuous limit: 

(V 1 (T,A*)> = 21n^V (2) 

where T is the duration of the Brownian path, and At is a regularization factor 
representing the minimum "Brownian-path time" along edges for them to be counted 
(see section CD). 

We then derive our main result, an exact formula for the average number of edges on 
the boundary of the convex hull of n independent Brownian paths of common duration T, 
starting at the same point (section 

(V n (T, At)) = a n (T, At) + n (T, At), (3) 

where a n (T,At) and (3 n (T,At) can be expressed as triple integrals. 

Finally, as a confirmation and illustration of our main formula (|3]), we show how the 
underlying probability distribution allows one to retrieve the existing formula for the 
average perimeter of the global convex hull of n independent Brownian paths (section EJ). 



1. The n — 1 - case 



For T > and t £ [0, T], let B(t) = (u(t),v(t)) denote the position of a two-dimensional 
Brownian motion at time t, that is: the coordinates u(t) and v (t) of B in the canonical 
basis of M 2 are independent, standard one-dimensional Brownian motions. They satisfy 
the simple Langevin equations: 

du , , dv , . / 1\ 

- = „„( f ), - = „„ (i ), (4) 

where r] u ,r] v are independent copies of Gaussian white noise with null average and 8- 
correlation, ie for all t, t! £ [0, T]: 

(Vu(t)} = (v v (t)) = o 

(Vu(t) Vu (t')) = (Vv(t)ri v (t')) = 5(t - t') , (5) 
(Vu(t)Vv(t')) = o 

S being Dirac's delta function, and (. . .) standing for the average over all realizations. 
We set B(0) = (0,0) and write Tb = B([0,T]) for the Brownian path and dCi for the 
boundary of its convex hull. 

Baxter's derivation of formula (CQ) relies on a combinatoric lemma that can be 
interpreted as a computation of the probability that the line segment joining two points 
of the walk separated by k steps belongs to the boundary of the walk's convex hull. 
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Figure 2: Given a Brownian path of duration T starting from the origin O, and given 
times s £ [At, T] and r £ [0,T — s], T B may be found to: (a) lie on one side only of 
the line through B(t) and B(t + s); or (b) cross it. In ), the line segment (solid 

black) joining B(t) and B(t + s) shall belong to the boundary of the path's convex hull 
; in case (b), it shall not. 

Similarly, in the case of a single Brownian path, we pick s and r, such that r £ [0, T] 
and (r + s) £ [0,T], and we wish to compute the probability for the line segment 
joining B{t) and B{r + s) to belong to dC\ (note that for s > 0, one has almost surely 
B(t) 7^ .B(t + s)). Among all possible realizations of a planar Brownian path starting 
at the origin and having duration T, some will indeed be such that the line segment 
[B(t)B(t + s)] appears on the boundary of the path's convex hull. Let us call such 
realizations (s, r)-paths. Using the Wiener measure for planar Brownian paths, we wish 
to compute the measure of the subset of (s, r)-paths, that is, the proportion of such 
paths among all possible realizations. Before we start this computation proper, let 
us observe that s represents the "Brownian-path time" between the two endpoints of 
the line segment under consideration. This time interval can be used to regularize the 
number of edges by requiring that: 



where At is some small, positive constant allowing one to actually "cut off' in our 
counting those segments that join points reached within too small a time lapse by the 
Brownian particle. 

Thus, once the parameter At has been chosen, our question becomes: for fixed 
s £ [At, T] and r £ [0,T — s], what is the probability that the whole path Tb lies 
on one side only of the line through B(t) and B{r + s), ie what is the probability of 
having figure [2a] rather than figure [2bl ? 

Let us examine r# in the coordinate system shown in figure [3aJ the line through 
B(t) and B(t + s) defines the x-axis, and the line perpendicular to it and through B(0) 
defines the y-axis. Note that if r > 0, B(0) = yo ^ almost surely, so one can fix the 



s £ [At,T], 
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Figure 3: The (s, r)-adapted coordinate system (a): the line through B{t) and B{r + s) 
defines the x-axis, and the line perpendicular to it and through B(0) defines the y-axis. 
If [B(t)B(t + s)] belongs to the boundary of the path's convex hull, the ^/-coordinate 
will be a combination of Brownian meanders and excursions (b). 

orientation of the axes by choosing -B(O) to indicate the positive or negative side of the y- 
axis. By property of isotropic planar Brownian motion, the new coordinates x(t) and y(t) 
of the path remain independent one-dimensional Brownian motions. We shall call this 
coordinate system the (s, r)-adapted coordinate system, because of the correspondence 
between planar Brownian paths for which the line segment [B(t)B(t + s)] belongs to 
the boundary of the path's convex hull, and pairs of linear, standard Brownian motions 
(x,y), where x is unconstrained except for x(0) = 0, while y satisfies the following 
constraints, as illustrated in figure l3bl 

CI. y starts at yo > (say) and hits for the first time at t — r; 

C2. y hits again at t = t + s but remains strictly positive between t — r and t = r + s; 
C3. y does not hit after time t = r + s, it remains strictly positive up to time t = T. 

Now, reasoning in the other direction, one can obtain a planar Brownian path Tb 
for which the line segment [B(t)B(t + s)] belongs to the boundary of the path's convex 
hull by picking two independent linear Brownian motions x and y behaving as stated 
above. Since the probability measure for planar paths is just the product measure 
derived from the measure for linear paths, one can compute the probability content of 
the subset of (s, r)-paths as the product of the probability contents of the corresponding 
x and y motions. 

The factor corresponding to the x-part is simple: any x-motion satisfying x(0) = 
is acceptable, so this will contribute a factor of 1 in the joint probability measure. 

As for the factor corresponding to the y-motion, it can be decomposed as the 
product of three terms thanks to the Markovian nature of Brownian motion (s and r 
being fixed): 

(i) the probability that the first part of the ^/-motion has a first-passage time equal 
to r (or equivalently, in reverse time, that it is a Brownian meander of duration r); 
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(ii) the probability that the second part is a Brownian excursion of duration s; 

(iii) the probability that the third part is a Brownian meander of duration T — (r + s). 

(Meander and excursion refer to the usual terminology for various types of constrained 
one-dimensional Brownian motion [441446] .) 

The first and third of these terms are readily found in the literature, since both can 
be viewed as the probability that a Brownian motion starting at some positive point 
first hits at a given time [4"5ll4"TH4"9] . Following [15], one can also write the second 
term from first-passage time densities. Indeed, if one introduces an intermediate time 

[t, t + s], the probability that y is an excursion between r and r + s can be split 
at k: after time k, y should first hit at time r + s; before time k, y should have last 
hit at time r (this is a last-exit time, as in [50l[5T]. Conditioned on y(n) = r > 0, 
this is equivalent to picking two Brownian motions starting at r, requiring that they 
have first-passage times n — r and r + s — k, and then putting them together back 
to back. Conveniently, the intermediate time k vanishes when one integrates over the 
intermediate position r = y(n). 

The threefold probability for the y-motion should be carefully normalized by that 
of a standard one-dimensional Brownian motion z subject not to the same constraints 
of course but to the same conditioning as y: z should be "pinned" at the same points 
z{t) = z{t + s) = 0, z(0) should be positive, and, most important, z(k) should be 
positive too. 

We are now in a position to write f(s, r), the measure of the subset of (s, r)-paths, 
in terms of p(t; u, v) and g(t; 0, v), where: 

p(t;u,v) = ^^exp [-( v -u) 2 /2t] (6) 

is the standard Wiener measure for a Brownian motion starting at u and reaching v at 
time t (this density is seen physically as the propagator of a free Brownian particle on 
a line); and: 

g(t;0,v) = -^=exp[-v 2 /2t] (7) 

is the measure of a Brownian motion starting at and reaching v > at time t without 
re-hitting in-between (see eg jl5]: this is the same as the first-passage time density 
for a Brownian motion starting at u > 0, and [g(t; 0, u) / / °° g(t; 0, v) dv] is indeed the 
propagator for the Brownian meander). 

Inserting the normalization detailed above and integrating over the initial, 
intermediate and final positions, one obtains: 

,/ x_ 9 ip 00 g(r;0,yo) dy J °° g(K-T;0,r)g(T + s-K;0,r)dr 



Jo 00 P( T > z o, 0) dz J °° p(n - r; 0, r') p(r + s - k; r', 0) dr' 
Jo°° g( T -(t + s);Q, y T ) dy T 
!Zv( T -{r + s)-^ZT)dz T 



(8) 
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Note that we have emphasized in the writing of the normalization the fact that it 
corresponds to a motion subject to the same conditioning as y: hence the different 
ranges of integration in the three parts of the normalization. Note also that we set 
arbitrarily y > in the constraints on the ?/-motion, but there are two sides to the line 
through B(t) and B(t + s), contributing an overall factor of 2 when moving from the 
probability measure of linear realizations to that of planar realizations. 

Incidentally, the middle term in (JSJ) can be interpreted as the probability that a 
planar Brownian path of duration s lies on one side only of the line joining its endpoints: 
this is 2ds/s, which is to be compared to 2/k, the probability that a discrete-time 
random walk formed with k > 2 steps lies on one side only of the line joining its initial 
and final positions |^2] . 

Now, performing the integrations in (jSJ) yields: 
iis^t[{T - s) - r\ 



And finally, f(s, r)ds dr represents the probability for [B{t)B{t + s)} to belong to dC\. 
Therefore, summing over all possible line segments with s > At: 

(V 1 (T,At)) = [ [ f(s,r)drds 



At JO 
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(10) 

as expected by taking the continuous limit of Baxter's formula (TjQ). 

Before we proceed to the general case, let us remark that (GO) can be interpreted as 
the product of: (i) a term corresponding to the probability that a Brownian motion of 
duration s is an excursion and (ii) a time-of-maximum density. This can be understood 
as follows: if one discards the excursion of duration s between r and r + s and joins 
the two parts left, one obtains a Brownian motion of duration T — s constrained to 
remain positive except at time r where it hits 0. Conditioned on its existence, the 
distribution of this time r is the same as that of the time at which a Brownian motion 
of the same duration attains its maximum — which is given by Levy's famous arcsine 
law [I], f(r) = 1/ ^y/r[(T-s)-T] 

2. The general case 

We now consider the global convex hull of n > 2 independent planar Brownian paths 
Bi, B2, . . . , B n , all starting at the same point and having the same duration T, as in 
figure [U Let us first observe that edges of the global convex hull will be of two types: 
edges joining two points of the same path r#- (type a) and edges joining two points 
belonging to two distinct paths T^., Tb with % ^ j (type b), as shown in figure HI 
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Figure 4: A sample realization of the convex hull of n = 9 Brownian paths in the plane, 
as in figure [TJ Edges on the boundary (dotted black line) of the convex hull can be 
of two types: (a) joining two points belonging to one of the paths (here for instance 
two points on B 9 (solid black line in the bottom right-hand corner)), or (b) joining two 
points belonging to two distinct paths (here on B 3 and B 7 (solid black line in the top 
left-hand corner)). 

Therefore, one has to consider not only line segments of the form [Sj(r)i?j(r + s)] with 
i G {1,2, ... ,n}, s G [At, T] and r G [0,T — s], but also line segments of the form 



[BiirjBjfa)] withz^ j G {1,2,..., n}, n G [0,T], r 2 G [0,T] and (n+r 2 ) G [At, 2T}. 



The reasoning is very similar to that in the n — 1 - case: 

(i) for type-a edges, we work in the (i, s, r)-adapted coordinate system, which is simply 
the (s, r)-adapted system for the i th motion. The ordinates of the n paths are linear 
Brownian motions y\, y 2 , ■ ■ ■ , y n which all start at the same point yo > 0; yi satisfies 
the same constraints C1-C3 as y in section [TJ and the (n — 1) other Brownian 
motions are simply constrained to remain strictly positive up to time T as shown 
in figure EaJ 

(ii) for type-6 edges, we work in an u (i,j, T\, T2)-adapted" coordinate system (the x-axis 
of which is defined by the line through Biiji) and Bj(r 2 )). yi,y 2 , ■ ■ ■ ,y n start at 
yo > 0; yi hits only once at time t — T\ and yj only once at time r 2 , while the 
in — 2) other motions remain strictly positive throughout, as shown in figure l5bl 

Again, the idea is to compute the probability of the 2-dimensional configurations 
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Figure 5: The ^/-coordinates of the n = 9 paths shown on figure [T] in the coordinate 
system adapted (a) to the edge of type a shown on figure 0] and (b) to the edge of type b 
shown on figure HI In both cases, all n motions start from the same positive value. In 
the system adapted to a type-a edge, (n — 1) motions are simply constrained to remain 
positive, while the one motion corresponding to the edge satisfies the same constraints 
as in section [I], shown in figure [3b] it hits at times r and r + s, and remains positive 
otherwise. In the system adapted to a type-6 edge, the 2 motions corresponding to the 
edge hit just once each (at times T\ and r 2 ) and remain otherwise positive, while the 
other in — 2) motions simply remain positive. 

by "counting" allowed one- dimensional configurations, for fixed s and r in the case of 
type-a edges, and for fixed T\ and r 2 in the case of type-6 edges. 

Using independence of the n motions, one can express the corresponding 
probabilities in terms of known propagators: p and g as defined in section [TJ along 
with q(t;x,y) = p(t;x,y) — p{t; —x,y). (q is the Brownian propagator between x > 
and y > in time t under the constraint of not crossing 0). One obtains: 

(i) for type-a edges (there are n paths to choose from): 



this is a term of a form very similar to ([8]), the only difference being in the first 
factor, which incorporates the probability that the (n — 1) other y-motions start 
from the same point y^ and remain positive throughout their common duration T; 

(ii) for type-6 edges (there are n(n — 1) ordered pairs of paths to choose from): 



ft\s,r) = 2xnx 



Jo°° g(T;0,y ) [/ °° g(T;y ,v)dv] n 1 dy 
jo°° p(t; zb, 0) p(T) z , v') dv' dz { 







x 



jlT 9(* - r ; 0» r ) 9(t + s- k; 0, r) dr g(T - (t + s);0,u) du 

fo° p( K ~ r i 0, r') p(r + s — k; r', 0) dr' p{T — (r + s); 0, u') du 




ti b \s, n, r 2 ) = 2 x n(n - 1) x 5(n + r 2 - s) x 
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Jo°° #( r i;0,?/oMT2;0,?/o) l/o°° <l( T ;yo,v)dv] n ~ 2 dy 



n ra-2 

dz 



o 

x 



/ °° p(n; zo, 0)p(r 2 ; z , 0) p(T; zb, u') dv' 
Jo 00 g(T-rt;0,u)du j °° g{T -r 2 ;0,w)dw 



Ho P(r - n; 0, u>) du> p(T - r 2 ; 0, w>) dw>'' 

this term is of course specific to n > 1 ; the ^-function in the first line simply ensures 
that the "Brownian-path time" along the line segment (one can go from end of the 
line segment to the other following a continuous Brownian path of duration n + r 2 ) 
is equal to a given s G [At, 2T] ; the following factors correpond to the probability 
that, having started at some yo > 0, Bi (resp. Bj) hits at T\ (resp. r 2 ) and 
remains positive before and after, while all other (n — 2) ^/-motions start also from 
yo and never hit for their entire duration T. 

Substituting the expressions for p, q and g in equations (TTTT) and (fl2"l) . tidying up 
the notations and changing to dimensionless variables in the integrals yield: 

pT pT—s pi pl—v poo 

a n {T,At) = / fl l a) drds= / / h^\u, v, z) du dzdv, 
J At Jo Jo Jo 



pT pT—s pi /"min(l,-u) poo 

/3 n (T,At) = ff ) drds= / h%\u,v,z)dudzdv, 

J At Jo J §± Jvaax.(0,v-1) J 



(13) 



(14) 



with: 



lrc-1 



wexp (-Jf) [erf(w)] ? 
/ii a) K^,z) = 4n ^ J - — (15) 

TTVZ [z(l — V — z)] 2 



and 



wVexp (--F^t) [erf(w)] n - 2 

h%\u,v,z) =4n{n-l) V — " } J r (16) 

nz{v — z) [irz(l — z)(v — z)(l — v + z)] 2 

where erf (u) = ^= J Q U e~ x2 dx is the error function. 

Finally, combining equations ( {TBI and (1T4"1) leads to our main result, an exact 
formula for the average number of edges ("cut off" at At) on dC n : 

(V n (T, At)} = a n (T, At) + n (T, At) (17) 

A first, simple check of this formula is readily carried out by looking at the case 
n — 2, for which the integrals over u are easily performed. We find: 

(V 2 (T, At)} = a 2 (T, At) + 2 (T, At) 

l)(l — v- z)\ 2 dzdv 

2 /■min(l,w) 1 

'max(0,i>— 1) 



7T I At 



2 />2 />min(l,i') 1 

+ -/ / [v 2 (l - z)(l - v + z)]~* dzdv 

7T J At 7 max Co,u-l) 
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= 2. n (g), (18) 

as expected indeed, since the case n = 2 is essentially the same, as far as the convex 
hull is concerned, as that of a single path of double duration. 



3. A new derivation of the average-perimeter formula 

Another confirmation of equation (ITT]) can be obtained by using the underlying 
probability distribution to retrieve the exact formula for the average perimeter (L n (T)) 
of dC n which was derived using an entirely different method in 



POO 

(L n (T)) = AnV^f / u e~ u * [eri(u)] n ~ l du. 
Jo 



(19) 



Indeed, multiplying the probability that a given line segment belongs to dC n by 
the average Euclidean length of this line segment before summing over all possible line 
segments should allow one to take the limit At = and re-obtain (TT9"j) . 

Let us therefore first compute the expected Euclidean length l(s) of the line segment 
[B(t)B(t + s)] (by the Markovian nature of Brownian motion, this average length is 
independent of r). This should not be done in the (s, r)-adapted frame, as this frame 
itself varies from one realization to another. Instead, one should use some fixed frame, 
eg the canonical basis of R 2 , in which the question now simply amounts to computing 
the average norm of a bivariate normal distribution with zero-correlation. This is a 
standard exercise and one finds that: 

l(s) = y^. (20) 

Note that for line segments of type b which are such that T\ + T2 = s, the average length 
also depends only on s and has the same expression. 

We shall now compute the expected perimeter length of dC n according to the 
probabilities established in section |2) 



(Ln(T)) 




1 rl—v r2 rm\n(l,v) 

l(vT)h^\u,v,z)dzdv + / l(vT)h%\u,v,z)dzdv 1 du 



max(0,u— 1) 



(21) 

Once multiplied by l{vT), the integrands in ( |2T|) become integrable from v = 0, meaning 
that we can set At = as expected. A few changes of variables then lead to a form 
where the integrals on v and z can be performed (see also [52]). One obtains: 

(L n (T)) = 4nV2rfT x 



POO POO 

/ ([erf^)]™" 1 - [erf( M )] n ) du+(n-l) / u (1 - erf(u)) e^Vf^)]™" 2 du 
Jo Jo 



(22) 
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A simple integration by parts of the second term in (I2"2"j) then gives: 

POO 

(L n (T)) = AnV^f / u e~ u ' [erf (u)]"" 1 du, (23) 
Jo 

exactly as in ( II 9p . thus confirming our main result, . 
Conclusion 

More work is needed to give (flTI) a more explicit form, but to the best of our knowledge 
this is the first exact formula for the number of edges on the boundary of the convex hull 
of planar Brownian paths. Progress on the treatment of the triple integrals appearing 
in ffl3|) and (I14p . including numerical evaluation and asymptotical approximations will 
allow comparison with results of computer simulations (which we have not included 
here) . 

Equivalent formulae for other types of stochastic processes, be they continuous-time 
or dicrete-time random walks in the plane, might be obtained by transposing some of 
the steps followed here. Further mathematical formalization of the probabilistic objects 
defined and used in this paper will be pursued and should lead to insights regarding 
potential extensions. 

Future work will also include looking at how the ideas introduced here, in particular 
the Brownian-path time "cut off", translate into the path-integral formalism which 
has proven quite powerful in contexts of constrained Brownian motions [53l[5l]- Such 
translation could be of interest especially in relationship with the study of Brownian 
cone points [6|[5T |[55T457] . 
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